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r — , We carry out a detailed stability analysis of the superconducting vortex solutions 

in the Weinberg-Salam theory described in Nucl.Phys. B826 (2010) 174. These vor- 
tices are characterized by constant electric current I and electric charge density Jo, 
*-~" 1 , for I —7- they reduce to Z strings. We consider the generic held fluctuations around 



the vortex and apply the functional Jacobi criterion to detect the negative modes in 
the fluctuation operator spectrum. We find such modes and determine their disper- 



sion relation, they turn out to be of two different types, according to their spatial 



behavior. There are non-periodic in space negative modes, which can contribute to 



S \ the mstab,lity of mfimtely lo,Ig vortlces ' but they can be el,mmated by ,mposi,ls 

the periodic boundary conditions along the vortex. There are also periodic negative 
modes, but their wavelength is always larger than a certain minimal value, so that 
they cannot be accommodated by the short vortex segments. However, even for the 
latter there remains one negative mode responsible for the homogeneous expansion 
instability. This mode may probably be eliminated when the vortex segment is bent 
into a loop. This suggests that small vortex loops balanced against contraction by 
the centrifugal force could perhaps be stable. 
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I. INTRODUCTION 

Superconducting strings (vortices) were introduced 25 years ago by Witten in the context 
of a simple field theory model containing two complex scalars and two Abelian vectors 



|l|. They generalize the well known Abrikosov-Nielsen-Olesen (ANO) vortex |2| for a non- 
zero longitudinal current supported by the scalar condensate in the vortex core [3J. The 
current can be very large, but there is an upper bound for it, which is typical for the 
superconductivity models, since too large currents produce strong magnetic fields which 
destroy superconductivity. Witten's superconducting strings have been much studied |3j, 
[4J, mainly in the cosmological context [5], J6|, since they can be viewed as solutions of 
some Grand Unification Theory [l] that could perhaps be relevant at the early stages of the 
cosmological evolution. 

The idea that superconducting vortices could also exist in the Weinberg-Salam theory 
was suggested longtime ago, because this theory, similar to the Witten model, includes scalar 



and vector fields and admits the 'bare' vortices - Z strings 



a 



It was therefore conjectured 



that it could also have 'dressed Z strings' containing a W condensate in the core 8|, |9|. 

10|, the 



However, when a systematic search of such dressed Z strings gave negative result 

€_idea was abandoned for many years. Only very recently it was reconsidered again 
and it was found that the negative conclusion of Ref. [lOj can be circumvented, 



who 





because it does not actually forbid the superconducting vortices to exist. Such solutions 



have been explicitly constructed in Refs. 



11 



12], but some of their properties turn out to 



be quite different as compared to those for the Witten strings. In particular, their current 
can be arbitrarily large, since its increase, although quenching the Higgs condensate, does 
not destroy superconductivity, because the current is carried by the vector W bosons and 
not by scalars as in Witten's model. 

Superconducting electroweak vortices can be viewed as generalizations of Z strings for 
non-zero electric current and charge. They exist for any value of the Higgs boson mass and 
for any weak mixing angle 6 W . Their current 1% and electric charge density I transform as 
components of a spacelike vector I a = (I , J 3 ) under Lorentz boosts along the vortex. The 
charge can be boosted away by passing to the 'restframe' where the electric field vanishes, 
while the current never vanishes and can be defined in the Lorentz-invariant way as X = 



a//| — Iq. The current is supported by the condensate of charged W bosons trapped in 
the vortex, while outside the vortex the massive fields die away and there remains only the 
Biot-Savart magnetic field produced by the current. In the X — > limit the vortices reduce 
to Z strings. For X ^> 1 they show a large region of size ~ X where the Higgs field vanishes, 
and in the very center of this region there is a compact core of size ~ 1/X containing the 



W condensate. In this estimates the unit of X corresponds to ~ 10 9 Amperes, so that the 
current can typically be quite large. These vortices could perhaps have interesting physical 
applications, but it is important to clarify their stability properties, which is the subject of 
the present paper. 

The fact that Z strings are unstable [131 ]. 14], [15] does not necessarily mean that their 



superconducting generalizations should be unstable too. The analysis in the semilocal limit, 
for 9 W = 7r/2, shows that the current-carrying vortices do possess instabilities, but the 
corresponding negative modes are all inhomogeneous, with the wavelength always larger 



than a certain minimal value depending on the current 16]. A a result, all instabilities can 
be removed by imposing periodic boundary conditions with a sufficiently small period. In 
this respect the vortex instability is qualitatively similar to the hydrodynamical Plateau- 
Rayleigh instability of a water jet, or to the Gregory-Laflamme instability of black strings 



in the theory of gravity in higher dimensions (see [17] for recent reviews) . 

Below we carry out the stability analysis for the generic superconducting electroweak vor- 
tices, for 8 W < 7r/2. We consider the most general field perturbations around the vortex and 
look for negative modes in the spectrum of the fluctuation operator. Our main conclusions 
are as follows. For any values of current X and charge I the vortex possesses inhomogeneous 
negative modes which can be periodic or non-periodic in space. These instabilities tend, at 
least as long as the linear perturbation theory applies, to split the vortex into non-uniform 
fragments. However, imposing the periodic boundary conditions with a period L along the 
vortex will remove all non-periodic negative modes. In addition, if L is small enough, all 
inhomogeneous periodic negative modes will be removed as well, similarly to what one finds 
in the semilocal limit 16j. Although this suggests that the periodic vortex segments should 
be stable 12], they still possess the homogeneous perturbation mode, which is not removed 
by periodic boundary conditions, since it can be viewed as periodic with any period. It 
is therefore important to know whether this mode is negative or not. It is actually non- 
negative for any X if 8 W = 7l"/2, and for any 8 W if X = 0, in which cases the periodic vortex 
segments can be stable. However, the detailed analysis reveals that for generic # W ,X the 
homogeneous mode is negative, so that the vortices remain unstable even after imposing the 
periodic boundary conditions. The instability makes them grow thicker. 

At the same time, it is possible that this remaining instability can be removed if the vortex 
segment is bent and its ends are identified to make a loop, since the thickness of a loop with 



a fixed radius cannot grow indefinitely. It is therefore possible that loops made of vortex 
pieces and balanced against contraction by the centrifugal force arising from the momentum 
circulating along them could perhaps be stable. This conjecture can be considered as the 
'positive' outcome of our analysis. Of course, its verification requires serious efforts, since 
one has to explicitly construct spinning vortex loops and then study their stability. However, 
any possibility to have stable electroweak solitons can be very important. 

The rest of the paper is organized as follows. In Sec. II the electroweak field equations 
are introduced and their vortex solutions are described. Sec. Ill considers the generic vortex 
perturbations, separation of variables in the perturbation equations, gauge fixing, and re- 
duction to a multi-channel Schrodinger problem. Sec. IV describes the Jacobi criterion used 
to reveal the existence of negative modes in the perturbation operator spectrum, as well 
as the explicit construction of these modes. The limits of zero current and large current 
are considered, respectively, in Sec.V and Sec. VI. The electrically charged vortices are dis- 
cussed in Sec. VII, while Sec. VIII contains concluding remarks. The two Appendices list the 
complete equations for the background fields and for their perturbations. 

II. SUPERCONDUCTING ELECTROWEAK VORTICES 

The bosonic sector of the Weinberg-Salam theory is determined by the action density 

C = ~~T2 W^W""" - -L B, U B^ + (D^D^ - ^ ($+$ _ l) 2 . (2.1) 

Here the Higgs field $ tr = (<&i, $2) is in the fundamental representation of SU(2), its covari- 
ant derivative is D^ = [d^ — ~ B hi — | r a W") $ with r a being the Pauli matrices, while the 
field strengths are W£„ = d^l - d^ + e abc W^W^ and B^ = d^B v - d u B^. The two cou- 
pling constants are g = cos 9 W and g' = sin 8 W where the physical value of the weak mixing 
angle is sin 2 9 W = 0.23. All quantities in (12.11) are rendered dimensionless by rescaling, their 
dimensionfull analogues (written in boldface) being B M = ^qB^, W^ = $ W^, $ = ^o*^, 
the spacetime coordinates x M = x M /go$o- Here $0 is the Higgs field vacuum expectation 
value and go relates to the electron charge via e = gg'hcgo- 

The theory is invariant under the SU(2)xU(l) gauge transformations 

$ ->■ U$, W -)■ UWIT : + 2tVd fl V~ 1 dx^ , (2.2) 



with U = exp (| i? + | r a 9 a ) where $, 9 a are functions of x M and W = (B^ + r a W^)c?x At is 
the SU(2)xU(l) Lie-algebra valued gauge field. Varying the action with respect to the fields 
gives the field equations, 



d> i B^ = g' 2 -((D 1/ <f>y<f>-&D v <f>), 

D^W a ^ = g 2% - ((A,$)V$ - &r a D v $), 







0. 



(2.3) 
(2.4) 
(2.5) 



1)^$ + - ($ f $ - 1)$ 

with DJNlp = dJNlp + CofeWjW^. 

The perturbative mass spectrum of the theory contains the photon and the massive Z, 
W and Higgs bosons with masses, respectively, being m z = l/\/2, Tt w = g^z, m H = \ffim z 
(in units of e& Q /(gg')). The exact value of the parameter /3 is currently unknown, but it is 
constraint to belong to the interval 1.5 < < 3.5. Defining the electromagnetic and Z fields 
as [IS | 

(2.6) 



Ffiu =- — B^ v - — n a W^ u , Z^ 



9 9 

with n a = &T a &/(&§) the electromagnetic current density is 



B, u + n a W% 



J„ = d u F, 



Vft ■ 



(2.7) 



A straight vortex oriented along the x 3 axis can be described by splitting the spacetime 
coordinates x^ into two groups: x k = (x 1 , x 2 ) spanning the 2-planes orthogonal to the vortex, 
and x a = (x°, x 3 ) parameterizing the 'vortex worldsheet'. Introducing the worldsheet vectors 



(sinh(fe), cosh(6)) 



(cosh(fe), sinh(fe)) 



U a — CTZjq, 



with b, a being two parameters, one makes the stationary, cylindrically symmetric field 



ansatz 



12| 



W = u(p) a a dx a — v(p) dip + r 1 [u\{p) a a dx a — i'i(p) dtp] 

h(p) 



+ t 3 [u 3 (p) a a dx a - v 3 (p) d(p], 



$ 



(2.9) 



where /i, J2 G R and the polar coordinates are introduced, x l + ix 2 = pe l(p . In what follows 
we shall call b, a, respectively, the boost and twist parameters. This ansatz keeps its form 
under Lorentz boosts along the x 3 axis whose only effect is to shift the value of b. This 



parameter is thus purely cinematic - one can always pass to the 'restframe' where b = and 
the field configuration is purely magnetic. The ansatz (12. 9p also keeps its form under gauge 
transformations ( 12.21) generated by U=exp{ — |IY 2 } with constant T, whose effect is 

(/1+1Y2) -> e F (/i+«/ 2 ), (ui+i« 3 ) -> e- ir ( Ml +iM 3 ), (ih+ivs) -> e^K+Wa)- (2.10) 
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Figure 1: Profile functions for the vortex solution with I = 2.57, n = v = 1, /3 = 2, sin 2 # w = 0.23. 

With the parametrization (12.91) the field equations (I2.3l) - (l2.5p reduce to a system of 
ordinary differential equations (IA.1[) — (1A.9D for the eight functions u, u±, u 3 , v , V\, v 3 , / 1; f 2 
listed in the Appendix A. It is worth noting that the boost parameter b drops from these 
equations, but they explicitly depend on the twist parameter a . The boundary conditions 
for the equations are obtained by requiring the energy density to be finite and the fields 
to approach at large p the purely electromagnetic Biot-Savart solution associated with the 
infinitely long electric wire. The local analysis in the vicinity of p = 0, 00 then gives the 
following boun dary conditions for the field amplitudes for ^— p — > 00 (keeping only the 



leading terms) 12| 



a\ <— u — > c\ + Q In p , In — v <— v — > c 2 , 

<— «i — > — (ci + Q In p) sin 7 , <— Vi — > — c 2 sin 7 , 

1 ^— « 3 — > — (ci + Q In p) cos 7 , v 4— V3 — > —C2 cos 7 , 
0^A^cos|, q SZ<- f 2 -+sml, (2.11) 

where 01, C\, c 2 , Q, 7, (? are real while n, v are integers. These boundary conditions imply in 
fact that the vector fields (12.91) are singular at the symmetry axis, but this singularity can 
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be removed by the gauge transformation (12.21) with U = e i(n v l 2 )^ e %vl P T I 2 which renders all 
fields ^-dependent, 

W={u{p)+T (p u 1 (p)+T 3 u 3 {p)} a a dx a (2.12) 



+ {2n - v - v(p) - r v vi(p) + t 3 \u - v 3 (p)}} dip, $ 



D i(n— v)<p 



hip) 



where t v = r 1 cos(vLp) — r 2 sin(i/if). The logarithmically growing at large p terms in the solu- 
tions is the specialty of the Biot-Savart field, which is essentially the Coulombian potential in 
two dimensions. The local analysis also shows 12| that the fields approach their asymptotics 
(12. lip for p — > oo exponentially fast as e~ miip , e~ mzP , e~ m " yP where m a = yff^ + a 2 u(p) 2 is 
the W-boson mass 'dressed' by the interaction with the long-range Biot-Savart field. 

Numerically integrating Eqs. (lA.l|) - (IA.8|) with the boundary conditions (12. lip gives the 
global solutions in the interval p G [0, oo) 12[. These solutions can be viewed as field- 
theoretic realizations of electric wires, where the wire is represented by a regular distribution 
of massive non-linear fields in the vortex core, while in the far field zone the massive fields 
die away and everything reduces to the pure Biot-Savart field. 

The vortices have the winding number n > 1 and the 'polarization' index v — 1,2,... u max 
where ^ ma x(/3, #w> n ) ranges from n for W = ir/2 to 2n — 1 for 8 W = 0. In addition, they are 
characterized by the worldsheet current vector 

2nQa c 



d p F an d 2 x 



fia 



99' 



XS r 



(2.13) 



Here J = Xsinh(6) is the electric charge per unit vortex length and J 3 = Xcosh(fe) is the 
total electric current through the vortex cross section. It is convenient to use X, b instead of 
I a as the solution parameters. 

The typical profiles of the solutions are shown in FigJTJ When X — > these solutions 
reduce to Z strings, that is to the embedded ANO vortices. For I ^ the amplitudes 
u, Ui,Us grow with p and show the logarithmic tails at infinity. Solutions up to X ps 12 



12] , and since the dimensionless value X = 1 corresponds to c<fr 



were constructed in Ref. 
1.8 x 10 9 Amperes, the vortex current can in fact be quite large. In addition, it seems that 
there is no upper bound for possible values of X, at least in the classical theory. This can 
probably be related to the fact that the vortex current is carried by the vector W bosons 
demonstrating the anti-screening effect 19| - the W condensate sets up currents which tend 



to increase the magnetic field and not decrease it as in the conventional Meissner effect. 



When considered in the restframe, where b = Iq = 0, the vortex is purely magnetic and 
characterized by its current X and the magnetic and Z fluxes. Charged (boosted) vortices 
with Iq 7^ have in addition the electric field, momentum and angular momentum. 

III. PERTURBING THE VORTEX 

Let us consider small perturbations around the vortex configuration (W", B^, $), 

W^W^ + SW^, B^^B^ + SB^, $^$ + <5$. (3.14) 

Inserting this into the equations (|2.3|) - (|2.5p and linearizing with respect to 5W":, SB^, 5$ 
gives the perturbation equations 

D^S® - i (5B^ + 5W*r a ) £P$ + - (2|$| 2 - l) 6$ + -5$ f $ 2 

= - (d^B^ + T a V„SW a ") $ , (3. 15a) 

a' 2 
d (l & M SB u + — {$ f (SB U + 8W av r a ) $ + 2i (5&D"$ - (£>"$) 5$) } 

= d u (d^B» + % 2— (5$t$ - $ t (5$) J , (3.15b) 

T> lt W i 6W av + e abc 5W*W c ^ 

a 2 
+ — {$V a <5$ (cJS" + 8W cv 8 a c ) + 2z (<J$ t r a £)' / $ - (£>"$) r a <5$) } 

= £>" (£> M cW^ + — (<5$V$ - $V a <J$) J , (3.15c) 

with P M X a = <9 M X" + e abc W*X c . 

These equations are invariant under the infinitesimal gauge transformations, 

<$$-K$$ + ^(to?+<J0 a T°)$, SBp-tdBp + dpSti, 5W*^5W* + V^6 a . (3.16) 

To suppress the pure gauge modes, we impose the background gauge conditions, 

d^B" + — (5$+$ - $ f 5$) = , 

2VW^ + ^|- (6&T a $ - &T a 5$) = , (3.17) 



which eliminates the right hand sides in Eqs. (13. 15b|) . (13. 15cp . However, this still lives the 



residual gauge freedom generated by parameters which fulfill the ghost equations (n a being 
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defined after q» 



g' 2 



d^5d + — $ f $ (&? + n a 66 a ) = 0, 
V^56 a + — $ f $ (n a cW + 58 a ) = 0. (3.18) 



A. Generic perturbation separation of variables 

Since the background fields depend only on the radial coordinate p, we can Fourier- 
decompose the perturbations with respect to x a ,ip. Keeping in mind the action of the 
Lorentz boosts on the background solutions, we wish to keep track of their action on the 
perturbations too. We therefore introduce S = (cjSq, + nT, a )x a + rrup, which reduces in 
the restframe to ux° + nx 3 + rrup. Denoting 5W® = SB^ the generic perturbations can be 
decomposed as 

5$ a = ^J {[0 a (<v, K,m\p) + iip a (uj, K,m\p)] cosE 

u,K,m 

+ [ir a (u, k, m\p) + i Xa(w, k, m\p)\ sin E} , 
-6W£iP = J2 {Xi(oo,K,m\p)cosE + Y 1 a (uj,K,m\p)smE} , 

u,K,m 

-5W^= J^ {X^ou, K,m\p) cosE + Y«(ou, K,m\p) sinE} , 

u,K,m 

<*Wk° = Yl W(u,K,m,p)coaE + Y^(u,K,m,p)amE} , (3.19) 



where a = 1, 2 and k — 1, 2 while now a = 0, 1, 2, 3. The infinitesimal gauge transformations 
Eq. ()3.16p can be decomposed in the same way (with 6° = $) 

S9 a = Y {«"(w, k, m\p) cos E + 7 a (w, re, m|p) sin H} . (3.20) 

Inserting the decompositions (I3.19P into Eqs.f )3.15|) the variables x a ,ip decouple and one 
obtains, for given u, re, m, a system of 40 ordinary differential equations for the 40 radial 
functions a , . . . , Y£ in f)3.19p . These 40 equations split (if re G M) into 2 independent 
subsystem of 20 equations each. These subsystems are identical to each other upon the 
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Ipa 


<r^ 


Xa j 


Y 2 

1 2 


•H- 


x 2 , 


X k 


-H- 


-Y 2 



replacement 



X 2 a ^-F 2 a , ^^-n 2 - (3-21 

Here a = 1, 2 but a = 0, 1, 3 and k — 1, 3, 4 (where possible we shall not write explicitly the 
arguments (u;, k, m|p)). Such a splitting of the equations into two groups is the consequence 
of the fact that the background configurations ( 12.91) are real, and so that the real and 
imaginary parts of their perturbations should be independent. In Sec lVHI below we shall 
study the case of complex k, u, and then the 40 equations do not split into two subsystems 
any more, but for the time being k is real and we can restrict our analysis to the 20 equations. 
These are equations for the 20 radial amplitudes in the right hand sides of (I3.2ip . they 
factorize with cos S. 

These equations are rather long and we do not write them down explicitly. Not all of 
them are independent, since there are four identities relating them to each other. These 
are the linearized versions of the identities obtained by taking the divergences of the vector 
field equations (12. 3p and (12. 4p . their existence is the manifestation of the gauge invariance. 
The equations are also invariant under the action of the gauge transformations, which now 
assume the following explicit form: 



x°-> 


lf-^7 , 






x\ 


-»■ 


Xl-cuj 1 , 


Y^ 


y 2 ° + (7°y , 






Y} 


-)■ 


Y 2 l + (7 1 )', 


Xl-r 


X 3 ° + m7 , 






x\ 


-)■ 


X3 + m 7 1 + v^a 2 , 


x 4 °-> 


X 4 ° + /,7°, 






x\ 


-» 


X A + k 7 1 — cru 3 a 2 , 


Y 2 ^ 


Y 2 + oo a 2 , 






x\ 


-)• 


Xl - u 7 3 , 


X 2 -^ 


X 2 + (a 2 )', 






Yi 


->■ 


^ 2 3 + (7 3 )', 


y^ 


Y 2 -m a 2 + 


(^i7 3 - 


- W37 1 ) , 


xl 


->■ 


Xl + m 7 3 - wi« 2 , 


Y 2 ^ 


Y 2 - k a 2 + 


o-(m 3 7 1 


- «i7 3 ) , 


xt 


->■ 


-^4 + « 7 3 + (JMia 2 , 



0i -> 0i + 2 « 2 / 2 , Xi ^ Xi + 2 [(7° + 7 3 )/i + 7V2] , 

02 -> 02 - i a 2 /i , X2 ^ X2 + ^bVi + (7° - 7 3 )/ 2 ], (3.22) 

where the prime denotes differentiation with respect to p. 
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B. Gauge fixing 



The additional terms on the right in (13.221) are pure gauge modes. They automati- 
cally fulfill the perturbation equations for any gauge functions 7 = 7°(w, K, m\p), 7 1 , a 2 , 7 s 
(verification of this is a good consistency check). We need to impose gauge conditions to 
eliminate these non-physical solutions. For example, one can use the temporal gauge, which 



completely eliminates all gauge degrees of freedom [16(. However, the fluctuation operator 
becomes then rather complicated. We have therefore chosen to use the background gauge 
conditions (I3.17p . they lead to more easy to handle equations, although not eliminating 
completely all gauge modes. 

After separating the variables the background gauge conditions (13. 17f) reduce to four 
constraint equations 



-uY 2 ■ 



(d p + -\ Y 2 ° + j 2 Xl + kXI + g' 2 (f lXl + f 2X2 ) = , 
(d p + -\ Y 2 X + j 2 X\ + kX\ + g 2 (f 2X i + /1X2) + au 3 Y 2 - ^Y 2 = 
d p + i) X 2 - ^Y 2 - kY 2 + g 2 (u, - fM 



a ( Ul Xi - u 3 Xl) + - (v 3 Xl - Vl X*) = 



ujXf 



d +-) Yi 



in 



Vl 



-XI + kX 6 a + g 2 (f lX i - /2X2) - a Ul Y 2 + ^Y 2 = . (3.23) 



pj - p* - p< 

Imposing these, one discovers that the 20 radial equations split into two independent sub- 
systems as 4+16, since the four amplitudes in (I3.23f) which are proportional to uj decouple 
from the remaining 16 amplitudes. Let us call these four amplitudes temporal, they are 
governed by the equations 



( D x S T\ 

S D 2 U W 

U D 3 V 

\T W V D A J 



I XV 9' \ 

X\/g 

Y?/g 

V X ll9 J 



(3.24) 
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where 
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m 



D 1 = --d p {pd p ) + ^ + ^-^ + *-(fl + fl), 
P P a 

D 2 = ~d p (pd p ) + ^±^a\l + ^-u> + 9 ^ (ft + ft) , 

p p l 2 v ' 

Ds = ~-d p (pd p ) + m2 + V l + V h \u\ + ul) + ^-^ + 9 ^ (ft + ft) , 
P P *• 

D, = --O p (pd p ) + ^±^a 2 ul + k 2 - ^ + t- {ft + ft) , (3.25) 

P P ^ 

and the off-diagonal terms are 

S = gg'hh, T = gg'(fl - fl) , 

U = -2 (^ + OT „ 3 ) , V = -2 (=1 - WMl ) , W = - (^ + A lM3 ) • (3.26) 

A direct verification reveals that if one resolves the constraints (I3.23P with respect to the 
temporal amplitudes, then the temporal equations will be automatically fulfilled by virtue 
of the equations for the remaining 16 amplitudes. The latter are described by Eqs. ()3.30|) 
below. Every solution of the 16-channel problem f)3.30p therefore generates a solution of the 
temporal equations (13. 24ft . This can be understood by noting that the temporal equations 
coincide with the ghost equations. 

The ghost equations describe the residual gauge freedom left in the background gauge. 
They can be obtained by inserting the pure gauge modes in ( 13.221) into (13.231) . or equivalently 
injecting the mode decomposition (I3.20p into (13. 18ft . This gives four radial equations for the 
gauge parameters 7°,7 1 ,a 2 ,7 3 which coincide with the temporal equations (13.241) upon the 
replacement 

Xl O 7° , X\ O 7 1 , Yl <B- -a 2 , X\ <B- 7 3 . (3.27) 

Therefore, the temporal amplitudes are pure gauge modes. It follows that they can be 
constructed via resolving the constraints (I3.23P if only the corresponding solutions of ( I3.30p 
are also pure gauge. Resolving the constraints for a non-pure gauge solution of (I3.30J) should 
also give a solution of the temporal equations (I3.24p . but since it cannot then be pure gauge, 
it can only be trivial. This gives a simple recipe to distinguish between the physical and 
unphysical solutions of the 16-channel Schrodinger system (I3.30P : if a solution fulfills the 
constraints ( I3.23P with zero temporal amplitudes then it is non-trivial, otherwise it is pure 
gauge. 
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We have explicitly tested this recipe for the negative modes of the system (I3.30p . These 
modes are all physical, since the spectrum of the ghost operator is positive, and because 
they fulfill the constraints (I3T231) with X° = X{ = Y? = Xf = 0. 

Since the four temporal amplitudes vanish for the physical solutions, one can use the 
constraints (13. 23ft in order to algebraically express four other amplitudes (for example those 
proportional to k) in terms of the remaining 12 amplitides. The system (I3.30P then reduces to 
12 independent equations only, which coincide with the equations obtained in the temporal 
gauge. However, their structure turns out to be rather complicated, which is why we prefer 
to work with the 16-channel system (13.301) . 

C. Reduction to a Schrodinger problem 

Imposing the background gauge conditions decouples the 4 temporal/ghost amplitudes, 
while the equations for the remaining 16 amplitudes can be cast into a Schrodinger form 
after the following operations. We redefine the amplitudes as 

Y 2 ° = ^=(^(Z + + Z^+A + + A^j , Yi = -L(g(Z + + Z_)-g'(A + + A-)), 
X°s = 9'^ (j (Z + - Z-) + A + - A^j , X* = -j= (g {Z + - 2L) - g' (A + - A-)) , 
X 4 ° = g' (V-Z + AA , Xl = (gZ - g'A ) , 

Yl = (yy+ + WZ + WZ + WZ) , Xl = - (W$ +WI-WI- WZ) , 

xl = p - (w+ - wz + wz - wz) , Yg = P - (-W+ + wz + wz - wz) , 
A = ^= (w - + w+) , Yi = -L {wz - w+) , 

02 = Yg i h * - ^ ' X2 = h ^ + ^ ' (3 ' 28) 

Here the notation A, Z and W^ 1 reflect the fact that these amplitudes correspond to the pho- 
ton, Z and W bosons, respectively. The subscripts refer to their polarizations. Introducing 
the 16-component vector 

* tr = (Z ,Z + ,Z_,Ao,A + ,A-,W+,WX,WZ,WZ,WZ,WZ,ht,h7,h+,hz) , (3.29) 
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the equations assume the form 

--(ptf'V +U{n,m\p)q = u 2 V, (3.30) 

P 

where U is a 16 x 16 symmetric potential energy matrix depending on the background fields. 
Its explicit form is given in the Appendix B. These equations are invariant under u — > —u, 
k — > — k and m —)■ —m provided that 

Z (u, n,m\p) —)• —Z (—lj, —k, — m\p) , 2±(u, k, m\p) — > Z t (—oj, —k, — m\p) , 

Ao(tu, k, m\p) — > — Ao(— to, —K, —m\p) , A±(co, k, m\p) — > A T (— to, —AC, — m\p) , 

Wq (w, n, m\p) — > — W c f(— u>, —K, —m\p) , W±(ui, n,m\p) — » W^(— ui, — ac, — m\p) , 
ht (w, ac, m|p) -> 7£(-w, -ac, -m|p) . (3.31) 

D. Boundary conditions 

The small p behavior of the perturbations can be determined by solving Eqs. (l3.30p in 
power series. For each of the 16 equations we find two solutions, one of which is bounded 
for p — y while the other one is divergent. The bounded solutions are 

Z v = c z J m -^ + ..., A n = c A J m ~^ + ..., W± = c f p \»±(™-v)\ + ..., 

hf = cfp lnTml + . . . , hf = cf p l"-^ m l + . . . , (3.32) 

where ct, c A c™ , c^ a are 16 integration constants and the dots stand for subleading terms. 

We are interested in bound state type solutions for which \1/ — > as p — > oo. In order to 

work out their behavior at large p, it is convenient to temporarily pass to the gauge where 

^2(00) = 0. This is achieved by applying the global symmetry (I2.10p with T = —7, which 

corresponds to the gauge transformation f)2.2p with U= exp{|7}. The background fields 

then simplify and one finds at large p 

h z h A h hl 

Z v = -5-e-"*' + . . . , A, = -^-q-^p + . . . , h+ + K = %z-» zp + . . . , 

y/P VP \fP 

b w± b h2 b hl 

W± = J^ e ~!^w ± d P + . . . , hf = -^ e -^ w ± dp + . . . , ht-h7 = —e~ PHp + ... . (3.33) 

' Vp Vp Vp 

Here the effective mass terms 



p 2 A = K 2 -u 2 , p 2 z = p\ + m 2 z , p\ = p\ + ml , p 2 w± (p) = (au(p) ± nf - u 2 + 



(3.34) 
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are assumed to be positive and b^, b*, W , b v a are 16 integration constants while the dots 
stand for the subleading terms. One can now apply to the whole system (background + 
perturbations) the inverse gauge rotation with U= exp{ — |7J. The background then returns 
to the gauge where /^(oo) = sin ^ while the perturbations (I3.33P change as 

Z v ->• (d' 2 + 9 2 cos 7) Z v + 2gg' sin 2 | A, - -J=Vl# sin 7 - -J=W" sin 7 , 
-At, -> (fl- 2 + g' 2 cos 7) A? + 2gg' sin 2 -^ + —= W+ sin 7 + — W~ sin 7 , 



W+ -»■ W+ cos 2 ^ - W~ sin 2 ^ + ^=Z r , sin 7 - -^= ^ sin 7 , 



2 7 ,.._ .j7, J 7 ■ 9 

^ 2"^ sm 2 + 7! ^ sm7 "72 

2 7 ,.,+ • 2 7, ~ • 9 

' 2" w " 8m 2 + 7! ^ sm7 "V2 



W~ -> >V~ cos 2 ^ - W+ sin 2 ^ + -^=Z V sin 7 - -^ A? sin 7 , 



/i* — y hf cos — — hf sin — , hf — > hf cos — + hf sin — . (3.35) 

This gives the large p behavior of perturbations. At this point we have everything we need 
to solve the perturbation equations ( 13.301) . 



IV. STABILITY TEST 

Summarizing the above analysis, we have arrived at the eigenvalue problem ( 13.301) and 
now we wish to know whether it admits bound state solutions with u 2 < 0. If exist, 
such solutions would correspond to unstable modes of the background vortex. In order to 
detect them, one possibility is to directly integrate the 16 coupled second order differential 
equations (I3.30p . However, if one just wants to know if negative modes exist or not, it is not 
necessary to construct them explicitly. A simple method to reveal their existence is to use 
the Jacobi criterion 



20|, which essentially uses the fact that the ground state wave function 
does not oscillate while the excited states do. It follows that if the zero energy wave function 
oscillates then the ground state energy is negative. 

A. Jacobi criterion 

When applied to our problem the Jacobi method gives the following recipe. Let \l/ s (p) 
with s = 1, . . . , 16 be the 16 linearly independent, regular at the symmetry axis solutions of 
( I3.30p . Each of them is a 16-component vector, ^ s {p) — ^f (p); I = 1, - - - , 16. Let A(p) be 
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the determinant of the matrix ^l(p 
part of the spectrum. According to 
of negative modes. 



. If it vanishes somewhere, then there exists a negative 



2l| . the number of zeros of A(p) is equal to the number 



Calculating A(p) is a much easier task than solving the boundary value problem (I3.30J) . 
since this simply requires to integrate the equations starting from p = with the boundary 
conditions ( I3.32p . This should be done, in principle, for each pair of values k, m. In 16| this 
method was used to test stability in the semilocal limit, where 6 W = n/2, while the typical 
behavior of the Jacobi determinant A(p) for W < tt/2 is shown in Figs j2|3l 





ln(l + p) 



ln(l + p) 



Figure 2: The behavior of the Jacobi determinant A(p) for fluctuations around the n = v = 1, 
I = 0.87 vortex (/3 = 2, sin 2 6 W = 0.23) for different values of k for m = (left) and for m = 1 
(right). The behavior for m = 2 is qualitatively the same as for m = 1. 

The main observation is as follows: the fundamental vortex with n — v — 1 has one 
negative mode in the m = sector for every value of n from the interval 

N < K max (J). (4.36) 

This can be seen in Figj2] where A(p) passes through zero exactly once if k is small and 
never vanishes if k is large. In fact, the symmetry relations (13.3 ip imply that u 2 (— k, —m) = 
lu 2 (k, m), so that for m = one has lo 2 (—k) = uj 2 (k) and it is therefore sufficient to consider 
only the k > region. In the Z string limit, for 1=0, one finds 

Nnax(0) = 2cr(0) (4.37) 

(see Table I) and also u 2 (0) = 0, so that the k = mode is not negative. For 1^0 one has 

K max (X) > 2<r(Z), (4.38) 
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and in addition we find that the k = mode is negative for 6 W ^ 7r/2, 



u; 2 (0)<0, X^O, 



(4.39) 



while in the semilocal limit one has uj 2 (0) = for all values of X 16 1. 

It seems that for the n — v — 1 vortex there are no other instabilities. We have checked 
for different values of k that for m — 1,2 there are no negative modes (see Fig]2]), while 
further increasing m increases the centrifugal energy thus rendering the existence of bound 
states less probable. As a result, it seems that the n — v — 1 vortices are unstable only in 
the m = sector and are stable with respect to any other perturbations. 
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ln(l + p) 
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Figure 3: The Jacobi determinant for fluctuations around the n = v = 2 vortex with X = 0.87 
{j3 = 2, sin 2 9 W = 0.23) in the m = (left) and m = 2 (right) sectors for different values of k. 



We also considered vortices with higher winding numbers n and v and found that the 
axially symmetric sector remains unstable for all solutions we examined (this was checked 
up to n = 3). An example is shown in Figj3]for n = v = 2. This instability is qualitatively 
the same as for the n — v — 1 vortex, it exists for < \k\ < K max (X). However, solutions 
with n > 1 have additional instabilities in sectors with m > 1 which can be interpreted as 
splitting modes. For example, the n = 2 solutions are also unstable in the m = 2 sector 
(see FiglS]), which apparently corresponds to breaking of the n = 2 vortex into two n = 1 
vortices. Such splitting instabilities are less interesting for us, and in what follows we shall 
concentrate on the intrinsic instability of the fundamental n = 1 vortex. 
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B. Finding the eigenvalue 



Having detected the negative modes, we now wish to construct them explicitly. Such 
a construction is considerably more involved than applying the Jacobi criterion, since it 
requires to solve the boundary value problem for the 16 coupled equations ( 13.301) with the 
boundary conditions (I3.32p and ( I3.35p . Unfortunately, even for m = these equations do 



not simplify much. We solve them with the multiple shooting method [22j, which requires 
to match at a fitting point the values of the 16 functions and their 16 first derivatives. It is 
then important to have enough free parameters in our disposal, and in fact we have the 16 
integration constants in the local solutions ( 13.321) . then 16 other constants in ( 13.351) . and also 
the eigenvalue u 2 . As we consider a linear system, one constant can be fixed by the overall 
renormalization, so that there remain 32 parameters to fulfill the 32 matchings conditions. 
Resolving these conditions gives us the global solution \P(p) of Eqs. ( 13.30]) in the interval 
p G [0, oo) and also the eigenvalue u 2 . 





-15- 1 = 3.55-^ 



1 = 4.13 



1 = 1.60 



1 = 1.94 



1 = 4.8 $- 



1 = 2.30 



-1.5 -1 



1 1.5 



Figure 4: Dispersion relation lo 2 (k) for the m = bound state solutions of Eqs. (|3.30p for the 



n = v = 1 vortex (/? = 2, sin 2 



0.23) for X < X* = 2.57 (left) and for X > X* (right). 



As a result, we obtain the dispersion relation u 2 (k) shown in FigJU We see that there is 
a value ft max (Z) such that lu 2 (k) < for |«| < K max . For small currents the function u 2 (k) 
has a double-well shape, with two minima of equal depth at n = ±n m i n and a local negative 
maximum at k — 0. As the current increases, K m i n decreases, the value u 2 (0) approaches 
a> 2 (±K min ), and finally K min vanishes for I = Z* when all three extrema of u) 2 (k) merge into 
a global minimum. For X > Z* the function u 2 (k) shows only one global minimum at AC = 0. 
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Some numerical characteristics of u 2 (k) are presented in Table [I] 




Figure 5: Profiles of K max (I) (left) and cr(I) against q(I) (right) for the same vortex solution as in 
FigfJ 



The passage from the two-well to one-well structure of the dispersion relation suggests 
that the system undergoes some kind of phase transition at X = X*. This is corroborated by 
the profile of K m i n (X) (see Figj5]) reminding of a second order phase transition. The point 
X = X* is also distinguished by the fact that the background 'consensate parameter' q = /2(C)) 
attains its maximal value there (see Figj5]). When X grows further, q starts decreasing and 
tends to zero as X — > 00. For large currents the vortex shows in its central part an unbroken 
phase region where the Higgs field is driven to zero by the strong magnetic field 12J. This 



suggest that the point X = X* corresponds to the transition in which the unbroken phase 
just starts to appear in the vortex center. The plot a(q) shows a characteristic two-branch 
structure (see Figj5]) and the point X = X* corresponds to the bifurcation between the 
two branches. Although this suggest that the stability may change at this point, we know 
already that the number of instabilities remains actually the same, but the dispersion relation 
changes its shape. As discussed in Sec JVHI below, this should alter the generic instability 
pattern. 

Let us now consider the limiting cases where the vortex current is either small or large. 
This will help to understand the structure of curves in FigJH 
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Table I: Parameter values for the n = v = 1 vortices with j3 = 2, sin 2 9 W = 0.23. 



X 


a 


c 2 (0) 


^min 


w ( K minj 


^max 





0.709697 


0.0 


0.709697 


-0.503670 


1.419394 


0.0804 


0.700 


-0.0370976 


0.705 


-0.519740 


1.415 


0.4851 


0.650 


-0.157024 


0.680 


-0.497821 


1.395 


0.8739 


0.600 


-0.255942 


0.655 


-0.502902 


1.390 


1.2430 


0.550 


-0.354806 


0.615 


-0.520995 


1.400 


1.6002 


0.500 


-0.462475 


0.570 


-0.558202 


1.425 


1.9494 


0.450 


-0.587058 


0.475 


-0.625174 


1.475 


2.3004 


0.400 


-0.738065 


0.280 


-0.741152 


1.560 


2.6740 


0.350 


-0.928761 


0.0 


-0.928761 


1.695 


3.0831 


0.300 


-1.12311 


0.0 


-1.12311. 


1.855 


3.5594 


0.250 


-1.44332 


0.0 


-1.44332. 


2.135 


4.1327 


0.200 


-1.89531 


0.0 


-1.89531. 


2.550 


4.8335 


0.150 


-2.56766 


0.0 


-2.56766. 


3.150 



V. ZERO CURRENT LIMIT 



Q, 



When the vortex current tends to zero, the solutions reduce to Z strings |7fl, whose 
stability has been studied before [13|, |l4|. The most detailed consideration of the problem 



was presented in Ref. 



13j, whose results we have been able to confirm. 



In zero current limit the vortex field amplitudes become 



u = —1, v = 2g' 2 (v ANO — n) + 2n — v = v z , u± — 0, u 3 = 1, 
vi = 0, v 3 = 2g 2 (v ANO - n) + v = v z3 , /i = / AN o = /z, h = °> 



(5.40) 



and the field equations (jA. lj) (IA.9[) reduce to the ANO system 



1 



(pfLo)' 



ANO 






1 ) I JANO 5 



Jano V 



ANO ""ANO 



(5.41) 
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whose solutions fulfill the boundary conditions <r- f ANO — > 1 and n <— v AN o -)■ as 
«— p — > oo. The solutions depend only on the winding number n, although when written 
in the gauge ( 12. 9 1) the fields also contain <r Q , ^, 




W z = (r d - 1) a a rfx Q - [v z (p) + T 6 v z3 (p)} dip, <$> z = | (5.42) 

The values of a 2 , v are determined by those for the generic vortices in the I — > limit, 
one has for example a 2 = a 2 ((3,6 w ,n,v) > 0. Although a a , v can be gauged away for this 
solution, they reappear again in the perturbation equations. In particular, —a 2 determines 
(see (I5.46P ) the eigenvalue in Eqs. (l5.43p . and since it is is negative, Z strings are unstable. 
Stable Z strings also exist, for unphysical values of (3,6 W (the eigenvalue is then positive), 



but they cannot be viewed as limits of superconducting vortices [12], so that they are not 
relevant for us. One can accurately determine the parameter regions in the (3, 9 W plane 
where Z strings are unstable/stable and so can/cannot be promoted to the superconducting 



vortices by studying solutions of Eqs. (l5.43p with a 2 = 13J, 12 1. 

Imposing (I5.40p . the potential energy matrix in the Schrodinger operator (IB.ip becomes 
block diagonal, so that the space of perturbations spanned by the 16-component vector \1/ in 
( IB.2P decomposes into a direct sum of six one-dimensional subspaces, one four-dimensional 
subspace, and two three-dimensional subspaces. The six one-dimensional subspaces are 
spanned by A±\, Aq, Z , W , which describe the photon and the longitudinal components 
of Z and W bosons. The potentials in the corresponding one-channel Schrodinger equations 
are positive definite so that there are no negative modes in these sectors. 

The four-dimensional subspace is spanned by Z±, h x , which correspond to the transverse 
components of Z and Higgs bosons. For m = this space further splits into sectors spanned, 
respectively, by Z + + Z-, h~l + h^ and by Z + — 2_, h^ —K[. Both of them contain bound 
states with u 2 > (in the first sector they exist only for (3 < 1.5) but there are no negative 
modes in this case. 

The remaining three-dimensional subspaces spanned by W+, h^ and VV_, h^ contain the 
negative modes. The perturbations are governed in this case by 

--{pV'±)' +U±^± = A±y± (5.43) 
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Figure 6: Dispersion relation (|5.47p for the bound state solutions of the eigenvalue problem (|5.43p 
(left). The two parabolas correspond to solutions in the independent $/ + and fy_ sectors. Passing 
to the gauge (J5.48P they get mapped into one parabola giving the dispersion relation for the modes 
(|5.49|) (right). The arrows indicate bifurcations with the superconducting branch. 



with A± = u 2 — (<t =F k) 2 and 



( 



u+ 



A™* V- 

A w * V T 



v ± y T ^ 



\ 



(5.44) 



where 



A w 



(2g 2 (v ANO -n) + u±(m- v))' 



A* 



[v ANO =f m 



2 (5 



± A V g 2V -^ + ^/ ANO , 



P 



<r 



+ a (Jano V + 9 ^ ANO ' 



V ± = g / ANO ± 



V ANOJ, 



ANO 



(5.45) 



p* 4 ' I \ p 

For m = equations (I5.43P admit bound state solutions both in the \l/ + and \l/_ subspaces 
with the eigenvalue 

A+ = A_ = -a 2 = -a 2 ((3, W , n, z/) < (5.46) 



for v = 1, . . . z/ max where rt < v max (/3, #w, n ) < 2n — 1 (l2[. These bound states are charac- 
terized, respectively, by the dispersion relation 



U! = UJ±{K) 



(a =f k,) 2 - a 2 = k(k =F 2a). 



(5.47) 



One has cj^(ac) < for < k < 2a and w£(«) < for — 2<r < k < so that there is 
one negative mode for every value of k from the interval (— 2a, 0) U (0,2a). As a result, 
the dispersion relation for negative modes is described by oj\{k) for k > and by w£(«) 



24 

for k < 0, therefore the uj 2 {k) curve consists of two parabolas intersecting at k = (see 
Figj6]). These parabolas continue to the \k\ > 2a regions where there are bound states 
with u 2 > 0. However, they should terminate for k = 0, since the effective photon mass 
p? A = k 2 — u 2 = ±2an defined by Eq. 03.34]) becomes imaginary after this point. Although 
the photon decouples for exactly vanishing background current, it rests coupled for however 
small but non-zero currents, when the background is arbitrarily close to Z string. 

Let us now use (I3.19P to reconstruct the dependence of negative modes on all spacetime 
coordinates. Then we apply to (I5.42p the gauge transformation U = e mif u(i'(p)u(a a x a ) with 
u(X) = e lX ( l ~ T )/ 2 . The Z string becomes then globally regular and independent of a a , i/, 

WT = 2(9' 2 + 9 2 r 3 )(n - v ANO (p)) *p, $T = ( e "^ o(p) J , (5.48) 

Q 

in which form it is usually described in the literature [Zfl. The negative modes read in this 
gauge (writing down only the Higgs field perturbations) <5<^> 2 = C + h^ {p)e {U)+ ^ k+ ^ t e~ %k+z for 
K e (0, 2a) and 5$ 2 = C_/i^(p)e^- (fc -^ t e* fc - z for k G (— 2a, 0). Here C± are integration 
constants, k± = k =f a and lu± = k\ — a 2 . Replacing k± — >■ k and using the fact that 
ht(p) = h^ip) one can wr he these solutions simply as 

5$ 2 = C±hj (p)e luj(k)lt e Tikz (5.49) 

with bJ 2 = k 2 — a 2 (see Figj6]). These negative modes can be viewed as standing waves of 
length A = 2ii/k whose amplitude grows in time. For k = ±er one obtains zero modes 
corresponding to the bifurcations of Z strings with the superconducting solutions. Since 
the minimal wavelength of negative modes is X min = 2ir/a, this suggests that the instability 
could be removed by imposing periodic boundary conditions along the z-axis with the period 
L < A m in - However, this would not remove the homogeneous k = mode, <5$2 = C^(p)e CTi , 
since it is independent of z and so can be considered as periodic with any period. 

Let us, however, consider Z string in yet another gauge - the one given by Eq. ( 12.121) . In 
this gauge the fields are also globally regular (we assume the restframe condition a a = a5^), 

W = a {r 3 - l)dz + W^ eg , $ = $ r z eg . (5.50) 

The correspondence between this gauge and (I5.48P is provided by the gauge transformation 
with 

r i 
U = u(az) = . (5.51) 

e laz 
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The negative modes ( I5.49P now become 

«J$ 2 = C±h+{p)e luJ±{K)lt e TiKZ (5.52) 

with k G (— 2cr, 0) U (0, 2a), these are standing waves of length A = 2tt/k > X min = ixjo. 
Imposing now periodic boundary conditions with period L = n/a will remove all negative 
modes. In particular, the mode which used to be homogeneous becomes now z-dependent 
with k = ±a, so that it will be removed. The ^-independent mode now corresponds to k = 
and it will not be removed, but this mode is not negative, so it is harmless. One should say 
that in the case under consideration all gauge invariant quanitites like 5B^ U and 5{n a W a ) 
vanish and there is no gauge invariant way to decide which modes are homogeneous. 

We notice finally that the gauge transformation ( 15.5 1[) is not periodic in the interval 
[0, vr/a], and therefore imposing the periodicity breaks the gauge equivalence between ( 15. 48ft 
and (I5.50p . The two descriptions of Z string become therefore physically different, which 
is why ( I5.50P becomes stable upon imposing the periodicity while (I5.48[) rests unstable. To 
the best of our knowledge, such a possibility to stabilize Z strings has never been discussed 
in the literature. 

Since the Z string zero modes for k = ±2cr correspond to bifurcations with the super- 
conducting solutions, they can be viewed as small deformations induced by the current. 
Now, one has u; 2 (±ft max ) = also for 1^0, which suggests that the related zero modes 
also correspond to deformations induced by a small current variation. However, for 1^0 
such deformations would inevitably contain logarithmically growing at infinity terms and 
therefore would not correspond to bound state solutions of the perturbation equations. This 
suggests that the k, = ±ft max zero modes could correspond to variations with respect to 
some other parameter. In other words, it may be that the vortex solutions admit stationary 
generalizations within a field ansatz more general than (12. 9ft . 

VI. LARGE CURRENT LIMIT 

When the current I is large, the vortex develops in its center a region of size ~ X where 
the magnetic field is so strong that it quenches the Higgs field to zero. Most of this region is 
filled with the massless electromagnetic and Z fields produced by the current. The latter is 
carried by the charged W boson condensate confined in the compact core of size ~ 1/X placed 
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in the very center of the symmetric phase. Outside the symmetric phase region the Higgs 
field relaxes to its vacuum value and everything reduces to the ordinary electromagnetic 
Biot-Savart field 12[. 



The vortex fields in this limit can be described by splitting the space into two parts: 
the core region p < Xq/X and the exterior region p > Xq/X. The fields in the core can be 
approximated by 

fx — h — <™3 = V\ = , uu = const. , v — 1 , 

a Ul (p) = XU X (Ip) , v 3 = V 3 (Xp) , (6.53) 

in which case the field equations (lA.lh — (1A.9[) reduce to 

lixUti^YLuu *(^y = ^V 3 , (6.54) 

with x = Xp. The solution of these equations exhibits the following behavior for <— x — > oo, 

0<-E/i(aO ^a\nx + b, 1 <- V 3 (x) -> (6.55) 

(here a = 0.29, b = —0.08 if g 2 = 0.23) where the large x asymptotic is attained, up to 
exponentially small terms, at x = Xq ~ 10. This determines the size of the core region. 
This solution describes the current-carrying charged W condensate confined in the core, the 
current value entering (16.531) as the scale parameter. 

The fields for p > Xq/X can be found separately and then matched to the core fields at 
P — Xq/X 12j. We do not need here the precise form of the p > x$/X solutions, since it is 
sufficient to analyze the stability of the core region. Indeed, suppose that we find a negative 
mode localized in the core. Since it vanishes in the outside region, it fulfills the perturbation 
equations also there, so that it will be a negative mode of the whole vortex configuration. In 
principle there could be additional negative modes in the outside region, however, it turns 
out that the core negative modes fit in well with the general instability pattern described 
above, which suggests that all vortex instabilities are localized in its core. 

To study the core instabilities, we inject (16.531) into the perturbation equations ( IB. II) . 
Passing to the radial variable x = Xp and defining 

Co = oo/X, k = k/X, (6.56) 
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the current X drops from the equations. For m = the equations split into three independent 
multichannel sectors plus free wave equations, and applying the Jacobi criterion one can 
check that the negative modes are contained only in the sector spanned by five amplitudes 



Y 2 

1 1 



X l (x), 



V2Xl = x{X z {x)-X A {x)) 



y3 

V2xl 



X 2 (x) , 

X 3 (x)+X 4 (x). 



X't = X 5 (x) , 



(6.57) 



Introducing the five-component vector and the potential energy matrix 
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x 2 
x 3 
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\x 5 J 
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( M x Q R \ 
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S T M_ U„ 
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(6.58) 



with the matrix elements 
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Q = -y/2d x Ux, R = -V2S = 
the unstable sector is described by 
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-2kU u U± = V2 [dJJx ± 



UxVz 



x 



M = Ul , 

U 2 



~{xV'X+UV = AV, 



(6.59) 



(6.60) 



where A = u 2 — k 2 . Fig|7] shows the Jacobi determinant A(p) for various values of k, and it 
seems that it always has a zero at some p, at least we could not find an upper bound K max 
beyond which A(p) ceases to vanish. Since such a bound always exists for small currents, it 
should presumably exist also for large currents, but to find it one should probably refine the 
approximation (I6.60P to take into account the region outside the core. At present, it seems 
that the description ( 16.601) is valid for any k if X — > oo or, if X is large but finite, up to some 
large but finite value of k. 

We then solve the eigenvalue problem ( I6.60P looking for bound states with the boundary 
conditions X 3 ~ X 5 = 0(1), X\ ~ X 2 = 0(x), X 4 = 0(x 2 ) at small x, while at large x 



X l ± X 5 ~ X 3 + X 4 T V2X 2 ~ exp{- f ^{U x =f k) 2 - Co 2 dx} 



X 3 — Xa 



exp{— v k 



Z-2 



uj 2 x) . 



.61) 
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This gives the dispersion relation lo 2 {k) shown in FigJTJ from where 



u\k)=1 i Cj\k/1). 



(6.62) 



We see that the negative mode eigenvalue is large for large currents, u ~ X, which relates 




ln(l + i) 




Figure 7: The Jacobi determinant for Eqs. (16.60j) for various values of k (left) and the dispersion 
relation uj 2 (k) (right) for the n = v = 1, /3 = 2, sin 2 6 W = 0.23 vortex in the large current limit. 



to the fact that the corresponding eigenmode is localized within a very short interval of size 
~ 1/X inside the core. It is worth noting that the one-well shape of this dispersion relation 
is qualitatively similar to what is shown in FigJH This suggests that the approximate 
description provided by Eqs. ( 16.601) is essentially correct. Since the ratio k/X is small unless 
K is very large, one has 

lo 2 (k) Ri X 2 u 2 (0) ps -0.12X 2 , (6.63) 

where the numerical coefficient is calculated for n = is = 1, f3 = 2, sin 2 9 W = 0.23. 

VII. CHARGED VORTICES 



Let us consider a bound state solution ^ K (p) of the Schrodinger problem (I3.30J) with 
the eigenvalue u 2 (k) (setting for simplicity m = 0). Injecting it into the mode decomposi- 
tion (I3.19P we reconstruct the dependence on all spacetime variables. The result will be a 
superposition of the real and imaginary parts of 



e iE ^ K (p). 



(7.64) 
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Here, using (j2.8p . 



S(£,z) = u b t + K b z 



(7.65) 



with (u;&, Kf,) being the Lorentz-transformed (boosted) components of the spacetime vector 
(w,«), 

Wf, = cosh(fr) a; + sinh(6) k, K b = cosh(6) k + sinh(6) uj. (7.66) 

The boost parameter is related to the electric charge J = Xsinh(6) (see (J2.13P ). Suppose 
that the mode under consideration is negative, u 2 < 0, so that u = i\u\. Then 



exp(i H) = exp {|o;|(cosh(6)t + sinh(&)z)} exp {in (sinh(6)t + cosh(b)z)} . 



(7.67) 



Let us consider first the uncharged vortex, for which one has b = I = and so 



exp(?S) = exp(|o;|t) exp(inz) 



(7.68) 



which grows in time but is periodic along z. Let us call such negative modes 'proper'. The ef- 
fect of this instability is schematically shown in Figj8]- the vortex undergoes inhomogeneous, 
periodic in z deformations which tend to segregate it into segments of length A = 2-k/k. Of 
course, this linear analysis is only valid as long as the perturbations are small, and so it does 
not imply that the vortex will actually break into segments - such a possibility is unlikely 
in view of the current conservation. Since the current density becomes inhomogeneous, this 
produces local inhomogeneities in the electric charge distribution in the form of a periodic 
sequence of positively and negatively charged regions along the vortex. 




/J = 3 




B = 2 



X=2nlK 



sin" ft, 



Figure 8: The effect of a proper negative mode on the vortex (left), and the eigenvalue lo 2 (0) of the 
homogeneous perturbation mode versus 9 W for fixed /2(C)) =0.1 (right). 
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The generic vortex perturbation can be decomposed into a sum over eigenmodes. As time 
increases, one can expect this sum to become dominated by negative modes whose growth 
rate is maximal (provided that the perturbation remains small). As is seen in FigJU |u;(k)| 
is maximal for k = ±ft m in if X < X* and for k = if X > X*. Therefore, for small currents 
the vortex will probably tend to segregate into segments of length 27r/ft m i n while for large 
currents it will rather expand homogeneously. 

Since the wavevector of all negative modes is bounded, \k\ < ft max , their wavelength is 
larger than 27r/fi; max . Therefore, imposing periodic boundary conditions with period 27r/K max 
will remove all these modes, because the vortex segment will not have enough room to 
accommodate them. Only the k = will stay, since it does not depend on z and so can be 
considered as periodic with any period. This poses no problems if X = 0, or if X ^ but 
6 W = it/2, since this mode is non-negative in these cases. However, for 6 W < n/2 and 1^0 
the homogeneous mode is negative. This is seen in FigJU i n Table I, and also in FigJH] which 
shows that o; 2 (0) is negative for X ^ and vanishes only for 8 W = n/2. This means that the 
generic vortices cannot be stabilized by periodic boundary conditions. 

We did not find any simple arguments explaining why w 2 (0) should generically be neg- 
ative. Since cj 2 (0) = for W = rr/2, when the massless fields decouples, one can suspect 
that the explanation could be related to the presence of the long-range field in the system. 
However, the massless fields decouple also for 6 W = 0, but in this case one has w 2 (0) < 
(see FigJH]). It is therefore likely that the explanation should rather be related to the non- 
Abelian nature of the background solutions. Indeed, the non-linear commutator terms are 
generically present in the backgrounds, but they vanish for 6 W = 7r/2 (when the SU(2) field 
decouples) or for X = (because Z strings are embedded Abelian solutions), that is exactly 
when cj 2 (0) vanishes. 

We have tried to analytically evaluate oo 2 (0) for small currents by using the method 
applied in Ref. J23J. In this method both the background and perturbation equations are 
expanded in powers of the small parameter q = /^(O) and then solved order by order. We 
have found that w 2 (0) = —cq 2 + . . . where c > 0, so that the homogeneous mode becomes 
negative for however small currents. It is also negative for large currents, as is shown by 
( I6.63p . Therefore, it is generically negative. 

Let us now consider electrically charged vortices with I ^ 0. Since they can be obtained 
by boosting the Iq = vortex, their perturbation can also be obtained in the same way. 
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Boosting the proper modes (J7.68P gives the negative modes (J7.67P of the charged vortex, and 
we shall call such modes 'boosted' in order to distinguish them from the proper modes. The 
boosted modes grow not only in time but also in space, along the vortex, which is a simple 
consequence of the fact that the time/space directions for the boosted vortex are not the 
same as for the Iq = vortex. The proper negative modes of the latter grow only in time 
when considered in the restframe, but the observer comoving with the charged vortex will 
see the very same modes grow not only in time but also in space (see Figj9]). Equivalently, 
one can say that the boost renders complex the wavevector k& in (I7.66p . since u is imaginary. 
Since the boosted negative modes grow with z, they can be used only within a finite 
range of z for the small perturbation theory to be valid. This can be achieved by forming 
wavepackets. Let us consider a wavepacket of the proper eigenmodes of the Iq = vortex, 



Sf(t,p,z) = / dKC{K)e^ t+iKZ V K (p) + 



(7.69) 



where the dots stand for the contribution of the scattering states (solutions of ( I3.30p which 
do not vanish at infinity). Assuming the initial perturbation 5f(0,p,z) to have a compact 




irolt 



- z 




Io =0 vortex 



Figure 9: Left: The proper negative modes grow along the spacetime history lines. The restframe 
time of the Iq = vortex flows in the same direction, while for the boosted vortex the time direction 
is different so that the instability grows not only with time t but also with z. Right: Spacetime 
evolution of the initial data with compact support S. 



support S along z-axis, its time evolution will be contained within the spacetime domain 
Y + (S) causally connected with S (see Figj9]). By simply changing the coordinates, t = 
cosh(b)t + sinh(b)z and z = cosh(b)z + smh(b)z, the same wave packet can be reexpressed 
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as sum over boosted modes, 

6f{t,p,z) = J dKC(K)e lE ^^ K (p) + ... (7.70) 

where H is given by (J7.65p ~ (j7.67p . If there are negative modes in (j7.69p . then (J7.70P will 
contain growing in z terms, but since z actually varies only within the finite range inside 
Y + (S) for a fixed t, the whole sum is bounded. One can therefore view this wavepacket as 
perturbation of the charged vortex with the initial distribution 5f(0,p,z) contained in S 
(see Figj9]). If Sf(t,p,z) grows with t then 5f(t,p,z) will grow with t, hence if the Iq = 
vortex is unstable then so is the Iq ^ vortex. 

So far, however, the symmetry between the I = and I ^ vortices is incomplete, 
because we have found only the proper negative modes for the former and only the boosted 
negative modes for the latter. These modes were obtained by solving the radial equations 
(I3.30p with real K and real u 2 < and they are spatially periodic in the vortex restframe 
but become non-periodic after the boost. One might therefore think that periodic boundary 
conditions could stabilize the charged vortices, since they will remove all boosted modes. 
However, there could be also solutions of Eqs. fl3.30p giving rise to negative modes which are 
initially non-periodic but become periodic after the boost. The boosted value «j in (I7.66P 
should then be real, hence one should look for bound state solutions of ( I3.30P for complex 
parameters 

o; = 7 -zfi, k = K + i£ltanh(b), (7.71) 

where 7 = j(b,K), Q = Q(b,K). It is worth noting that a similar recipe was consid- 
ered within the stability analysis of the boosted black strings in the theory of gravity |28|. 
Inserting this in (J7j66]), the imaginary part of Kf, vanishes and one obtains 

exp(i E) = exp (fi&£) exp (i^t + in h z) (7.72) 

with 

^■b = — tttt ' lb = cosh(6)7 + sinh(b)K, n h = cosh(b)K + sinh(6)7. (7.73) 

cosh(o) 

Since exp(z E) grows in time and has the harmonic z-dependence, this corresponds to proper 
negative modes of the charged vortex with Iq = Xsinh(fe). 

The next question is whether solutions of Eqs. (j3.30p for the complex parameter values 
(I7.7ip exist. If b = then k is real and we should recover the already known solutions with 
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0.25x11, 1 = 2.67- 




0.5 1 



Figure 10: Left: Real and imaginary frequency parts 7^ and fi;, versus Kb for perturbations of the 
boosted vortex with X = 0.08 for several values of the boost b S [0,0.4]. Right: 0, versus the 
boost parameter b for the homogeneous (k^ = 0) perturbation mode for the boosted vortices with 
X = 0.08 and X = 2.67. In both panels /3 = 2, sin 2 6 W = 0.23, n = v = 1, m = 0. 



real u 2 = — Q 2 , therefore one has 7 = in this case. If b 7^ then both u 2 and k become 
complex so that equations ( 13.301) should be complexified as well. We therefore obtain a 
system of 16 linear complex equations, which is equivalent to 32 real equations. This does 
not mean that the number of degrees of freedom doubles, since even for b = we actually 
had 32 (=40-8) real equations split into two independent subsystems of 16 equations each, 
equivalent to each other upon ( I3.2ip . For b 7^ we still have the 32 equations, but they no 
longer split into two subsystems. 

Solving numerically the 32 coupled equations is considerably more time consuming than 
solving the 16 equations. This is why we did not carry out a systematic analysis of the 
parameter space but studied instead just several representative cases. We integrated the 
equations looking for bound state solutions with the boundary conditions given by the 
complexified version of (I3.32p .f l3.33p . Choosing a value of b, we have managed to explicitly 
construct such solutions and determine 7 and Q as functions of K. 

It turns out that the dispersion relation for Qb against Kb remains qualitatively the same 
for b 7^ as for b = (see Fig JTQj) . If the current X is small and b is fixed then Q(Kb) 
starts at a non-zero value at Kb = 0, increases and reaches maximum, then decreases and 
vanishes for some maximal value Kb = ft m ax (see Fig flOl) . For large currents Q(Kb) decreases 
monotonously from its value at Kb = till zero. One has Q(Kb) = Q(—Kb). 
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The proper negative modes therefore exist for any value of charge and not only for I Q = 0. 
To completely restore the symmetry between solutions with different I , we note that the 
proper negative modes for any b can be boosted towards a different value of the boost 
parameter, B, say. This will give boosted negative modes of the Iq = Isinh(i?) vortex 
proportional to 



exp {Qj, cosh(5 — b)t + Qb sinh(£? — b) z} exp {i^Bt + k-bz)} 



(7.74) 



with Kb = cosh(B) K + smh(B)'-f and 7^ = cosh(i?)7 + sinh(i?)A'. Therefore, for any vortex 
charge J = Xsinh(B) there are proper negative modes, but also infinitely many boosted 
modes labeled by b ^ B. The space of negative modes has the same structure for any value 
of charge, since there is one-to-one correspondence between modes for different charges via 
boosts, as schematically shown in FigJTTl 




proper 
modes 



Figure 11: The set of negative modes for any given vortex charge, for example Iq = Xcosh(6), is 
represented by the vertical line. The proper modes are represented by the fat points. There is 
one-to-one correspondence between modes for different Iq via boosts. 



The boosted negative modes are non-periodic in space and can contribute only to the 
instability of infinitely long vortices, but they will be removed by imposing on vortex the 
periodic boundary conditions. The proper modes will stay then, but if the period is less 
than 2n/n max then they will be removed as well, apart from the Kb = mode. We know that 
for Iq = this mode is generically negative, but perhaps things may change for Jo 7^ ? 
We therefore trace Q for this mode against b and find that it decreases very rapidly with 
b, especially for small currents (see Fig JT0|) . so that the instability growth rate decreases 
when the vortex charge Iq increases. However, it is not clear from these data if Q always 
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stays finite or eventually vanishes at some large value of b. It seems however that the latter 
option is impossible, since Kb = implies that 7 = K = 0. Setting Q = would therefore 
mean that uj = k = 0, but since k = is real, this solution should be contained in the 
previously obtained dispersion relation u 2 (n). However, we know from the previous analysis 
that u 2 7^ for k = (unless X = 0), and so the value Q = is impossible. Therefore, there 
is no critical value of boost for which the homogeneous instability would disappear. 

VIII. CONCLUSIONS 

We study in this paper the stability of the superconducting vortex solutions in the 



Weinberg-Salam theory described in Ref. [12|. Such vortices are characterized by a constant 
electric current ^3 = Xcosh(fr) and linear electric charge density I = Xsinh(6) comprising 
a spacelike vector (Iq,I^). Fixing X, vortices with different values of the charge Iq can be 
related to each other by Lorentz boosts, in particular there exists the restframe where Iq = 0. 
For X — y all solutions become Z strings, while for 8 W — > tt/2 and /3 > 1 they reduce to the 
twisted semilocal strings studied in Ref. |24| |. 

We consider generic vortex perturbations in the linear approximation and find that af- 
ter separating the variables the perturbation equations reduce to the effective 16-channel 
Schrodinger problem f)3.30p . This problem admits bound state solutions with u 2 < whose 
dispersion relation u 2 (n) is shown in FigJH and tabulated in Table I. These solutions de- 
scribe the 'proper' negative modes of the I = vortex. Choosing the parameters u, K in 
Eqs. fl3.30p to be complex gives bound state solutions describing proper negative modes of 
the charged vortices. As a result, for any given value of charge Jo there is a one-parameter 
family of proper negative modes which can be labeled by the wavevector «&. These pertur- 
bation modes grow in time favoring segregation of the homogeneous vortex into segments, 
although one cannot conclude from the perturbative analysis whether it will actually break 
in pieces in the long run. 

Since vortices with different I are related by Lorentz boosts, their perturbations can be 
related in this way too. Boosting the proper negative modes of the I = Xsinh(6) vortex one 
obtains negative modes of the Jo = Xsinh(Z?) vortex, so that the latter acquires in fact an 
additional infinity of negative modes labeled by b ^ B. These 'boosted' modes grow with 
z but they can form localized wavepackets to contribute to the instability of infinitely long 
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vortices. Since they are non-periodic in space, they can be removed by imposing periodic 
boundary conditions along the vortex. 

The proper negative modes are proportional to exp{%Kbz} and so they can be made 
compatible with the periodicity along z by adjusting the value of Kb- However, they exist 
only for |/fy,| < K max and so choosing the period to be less than 2n/K raax the vortex segment 
will not have enough room to accommodate these modes. All of them will therefore be 
removed, apart from the Kb = mode which is independent of z and can be considered as 
periodic with any period. Therefore, the only remaining vortex instability is associated with 
this homogeneous mode. 

In some cases one has u = for Kb = 0, as for example for 6 W = n/2 and for any X 
(semilocal vortices), or for X = and for any 8 W (Z strings). In these cases the homogeneous 
mode is not negative and so the short periodic vortex segments turn out to be stable. In 
particular, Z strings can be stabilized in this way by passing to the gauge ( I5.50P and then 
imposing the periodic boundary conditions which break the gauge invariance. However, in 
the generic case the homogeneous mode is negative and it renders the vortex unstable with 
respect to the homogeneous expansion even after imposing periodic boundary conditions. 

At the same time, it is possible that the homogeneous negative mode could be removed 
by the curvature effects. Specifically, let us suppose that one 'cuts out' a finite vortex 
segment, bends it and identifies its extremities to make a loop. Then, since the loop thickness 
cannot be larger than its radius, the homogeneous expansion of the vortex segment should 
inevitably stop at some point. Therefore, the homogeneous instability will be removed, 
suggesting that the loop could be stable. Of course, this argument is only qualitative. 
Moreover, new instabilities could appear when bending the vortex. However, any possibility 
to have stable electroweak solitons, as for example vortex loops, could be very important. 
Such loops could be balanced against contraction by the centrifugal force arising from the 
momentum circulating along them. Since momentum flows along vortices with Iq ^ 0, they 
can be naturally used to 'make' the loops. All this suggests that spinning vortex loops - 
electroweak analogs of the 'cosmic vortons' (25| - could exist and could perhaps even be 
stable. Of course, verification of this conjecture requires serious efforts, since so far vortons 
have been explicitly constructed only in a simple scalar field model 26|, 27|. However, if 



the electroweak vortons indeed exist and are stable, they could be a dark matter candidate. 
There could be other physical manifestations of the superconducting vortices. They could 
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perhaps be created either at high temperatures or in high energy collisions, and since they 
are non-topological, they could exist in the form of finite segments. If their extremities are 
attached to something (charged clouds), then they could be spatially periodic and transfer 
charge between different regions of space like 'electroweak thunderbolts'. Non-periodic vor- 
tex segments should decay emitting jets of W ± through its extremities, which could perhaps 
be detectable at the LHC. Specifically, large magnetic fields similar to those inside the vortex 
and also large currents can be created in the LHC heavy ion collisions. This could lead to 
creation of virtual vortex segments whose subsequent disintegration would be accompanied 
by showers of W s. As a result, if one observes an excessive W ± production in the colli- 
sions, this could indicate the vortex segment creation. A similar way to detect the presence 
of the non-perturbative electroweak structures in the LHC collisions was discussed in |29[. 
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APPENDIX A. BACKGROUND FIELD EQUATIONS 

With the parametrization ( 12.9[) the field equations ( 12.3l) -( l2~5j) reduce to two U(l) equa- 
tions (with ' = j-) 

- p (puJ= ^-S( u + u 3 )f^ + 2uJJ 2 + (u-u 3 )A , (A.l) 



P[j)'= 9 -^U v + ^3)1? + 2vJJ 2 + (v -v 3 )f 2 \ , (A.2) 
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two Higgs equations 



- p (pfi)' = { y K«+ ^ + «?] + ^ [(« + ^3) 2 + Vl] + P-{fl + f 2 - 1)} f X 
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' y MMl + v Wl 



four Yang-Mills equations 
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P I — I = +<r 2 (ui«3 - V3U1) u 3 + j [vi(fi + A 2 ) + 2u/i/ 2 ] , 

f- J = -^ (Vl«3 " «3«l) Wl + y [(«3 + t>)/l + (V 3 ~ v)/f] 



and a first order constraint 
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'i u i u 3 - u z u 'i) + t ( u i v 3 - ^i) - # 2 (AA - AA) = 0. 
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(A.5) 
(A.6) 

(A.7) 
(A.8) 

(A.9) 



APPENDIX B. PERTURBATION EQUATIONS 



Fixing the gauge, decoupling the ghost modes as described in the main text and using the 
parametrization ( I3.28J) for perturbations in the physical sector, the perturbation equations 
can be written in the form of the 16-channel Schrodinger problem 
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where the 16-component vector \1/ and the 16 symmetric potential matrix U read 
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Here 
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The vector- vector couplings are defined by 
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where x and y design Z, A, W + , W and 

C = -^ (I + <^ 2 + (s 2 - g' 2 ) fl) , d 
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Finally, the vector-scalar couplings are 
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where 



a 



a l = 77 (( M + u 3)/l + M1/2) 



a 2 = 2 (( n ~ n3 )^ 2 + Ul ^ 



(B.10) 



,± = -L ^ ± i ((„ + <, 3 )/i + v 1/2)) , 4 = -^ U ± ^ ((« - t*)/ 2 + wx/O^ . 
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